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Abstract 



In this paper we preliminarily explore the possibility of designing 
a dedicated satellite-based mission to measure the general relativis- 
tic gravitomagnetic Lense-Thirring effect in the gravitational field of 
Mars. The focus is on the systematic error induced by the multipolar 
expansion of the areopotential and on possible strategies to reduce it. 
It turns out that the major sources of bias are the Mars'equatorial 
radius R and the even zonal harmonics Jg,£ = 2,4,6... of the areopo- 
tential. An optimal solution, in principle, consists of using two probes 
at high-altitudes (a w 9500 — 9600 km) and different inclinations (one 
probe should fly in a nearly polar orbit), and suitably combining their 
nodes in order to entirely cancel out the bias due to SR. The remaining 
uncancelled mismodelled terms due to 5J^,£ = 2,4,6, ... would induce 
a bias < 1%, according to the present-day MGS95J gravity model, 
over a wide range of admissible values of the inclinations. The Lense- 
Thirring out-of-plane shifts of the two probes would amount to about 
10 cm yr~L 

Keywords: Experimental tests of gravitational theories; Lunar, plane- 
tary, and deep-space probes; Mars; Gravitational fields 
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1 Introduction 

Recent claims concerning a detection of the general relativistic gravitomag- 
netic Lense-Thirring effect [1] in the gravitational field of Mars with the Mars 
Global Surveyor probe [21[3j|4l[5] raised interest concerning such a new Solar 
System scenario for testing Einsteinian gravity. In view of what happened to 
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the Gravity Probe B mission [61 [7] , which might finally reach an accuracjo 
far from the originally expected 1% in measuring the SchifF precessions [8] of 
the spins of four gyroscopes carried onboard, and of the lingering controversy 
about the realistic accuracy reached with the terrestrial LAGEOS-LAGEOS 
II test of the Lense-Thirring effect [5] (and reachable with the recently ap- 
proved LAGEOS-like LARES satellite, to be launched at the end o^| 2009 
by the Italian Space Agency [TU]) we feel that attempts to scour new and 
unexplored routes may be of some value. Indeed, as many laboratories and 
methods as possible should be used to extensively test a fundamental inter- 
action like gravity and its predictions, especially when their tests are so few 
and their outcomes uncertain. 

In this paper we wish to fix a first stick by studying in some details the 
possibility of using a dedicated mission to Mars with one or more probes to 
measure at a reasonable level of accuracy (a few percent) the elusive Lense- 
Thirring effect on the longitude of the ascending node 0, of the spacecraft's 
orbital plane. We will concentrate here on one of the major source of sys- 
tematic error, i.e. the multipolar expansion of the martian gravitational 
potential which induces on the node a huge noise with the same tempo- 
ral signature of the relativistic signal of interest (a linear rate) in order 
to see what are the critical issues in view of the present-day knowledge of 
the Martian space environment. We will not discuss here the perturbations 
of non-gravitational origin which depend on the shape, the instrumentation 
and the orbital maneuvers of such probes. They would be strongly related to 
possible other tasks, more consistent with planetology, which could be fruit- 
fully assigned to such a mission in order to enhance the possibility that it 
may become something more than a mere, although-hopefully-interesting, 
speculation; in this respect the medium-long term ambitious programs of 
NASA to Mars may turn out to be useful also for the purpose discussed 
here. 

1 See on the WE B |http: / / einstein.staJiford.edu/| 

2 See on the WEB http://www.esa.int/esapub/bulletin/bulletinl35/bull35Lbianchi.pdf 
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2 The use of one nearly polar spacecraft 

The Lense-Thirring effect consists of a small secular precession of the nod^l 
Q of the orbit of a satellite moving around a central slowly rotating mass 

■ _ m 

LT " C 2 a 3 (1 _ e 2)3/2' W 

where G is the Newtonian constant of gravitation, S is the spin angular 
momentum of the central body, a and e are the semimajor axis and the 
eccentricity, respectively, of the satellite orbit. According to the latest de- 
terminations of the global properties of Mars pXJ, S = (1.92 ± 0.01) x 10 32 
kg m 2 s _1 for the red planet. 

The oblateness of the central body, of mass M and equatorial radius R, 
makes the node to secularly precess as well according to [12] 

siw=^n d j t , (2) 

1=2 

where the coefficients explicitly computed in, e.g., p3j up to degree 
I = 20, depend on the orbital semimajor axis, the eccentricity and the 
inclination i to the planet's equator; for example, for £ = 2 we have 

• 3 f R\ 2 cosi 

where n = a/ GM/a 3 is the probe's Keplerian mean motion. Although 
the multipolar expansion of the classical part of the gravitational potential 
is usually included in the dynamical force models used to process probes' 
data, such a modelling is always imperfect; thus, an aliasing trend due to 
the mismodelling in the even zonals Jt and the other parameters like R and 
GM entering the classical node precessions may affect the recovery of the 
relativistic precession of interest. 

2.1 The even zonal harmonics 

Let us, now, focus on the systematic error 5fi induced by the uncertainty in 
the even zonals 5Ji by assuming for them the calibrated covariance sigmas 
of the latest Mars gravity model MGS95J [TTJ. We will evaluate it as 

Snj t <\n.t\SJt, £ = 2,4,... (4) 

3 Also the pericentre uj of a satellite secularly precesses under the action of the gravit- 
omagnetic force, but we will not consider here such an effect. 
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Figure 1: Systematic percent errors 5/j,j e per degree £ due to the mismod- 
elling in the even zonal harmonics SJi, £ = 2, 4, 6... for a = 3796 km, e = 0.01, 
89 deg < i < 91 deg according to the calibrated covariance sigmas of the 
MGS95J global gravity solution up to £ = 20. The Lense-Thirring effect 
amounts to 34 mas yr _1 corresponding to a shift in the out-of-plane direc- 
tion of 0.62 m yr~ x . 

Because of the presence of cos % in all the coefficients ti g of the oblateness- 
induced node precessions, we will concentrate on nearly polar (i « 90 deg) 
orbital configurations to minimize such a corrupting effect. 

It turns out that altitudes of a few hundreds of km typical of the ma- 
jority of the currently ongoing martian missions are definitely not suited 
for our scope, as shown by Figure [TJ Note that inserting a probe into an 
areocentric orbit is not an easy task, so that we decided to allow for a depar- 
ture of up to 1 deg from the ideal polar orbital configuration to account for 
unavoidable orbital injection errors. Also typical mission requirements pull 
the inclination some degrees apart from 90 deg: for Mars Global Surveyor 
(MGS) i = 92.86 deg. The systematic bias per degree increases for higher 
degrees £ = 2, 4, 6, ... and it turns out that only a very narrow range for i, 
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a= 9500 km, e = 0.01 
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Figure 2: Systematic percent errors 5/j,j e per degree i due to the mismod- 
elling in the even zonal harmonics SJi, £ = 2, 4, 6... for a = 9500 km, e = 0.01, 
89 deg < i < 91 deg according to the calibrated covariance sigmas of the 
MGS95J gravity model up to t = 20. The Lense-Thirring effect amounts to 
2 mas yr -1 corresponding to a shift in the out-of-plane direction of 0.10 m 
yr" 1 . 

i.e. Ai ~ 10~ 3 deg, unlikely to obtain, might push the systematic errors 
below the 1% level. 

By keeping a near polar geometry, much larger values of the semimajor 
axis, comparable with that of Phobos (a = 9380 km), one of the two natural 
satellites of Mars, yield reasonable results, as shown by Figure EJ Indeed, 
now the noise decreases with high degree terms and the most serious contri- 
bution is due to the first even zonal I = 2 yielding a bias up to about 11% 
for Ai = 1 deg; J4 affects the the Lense-Thirring precessions at most at 3% 
level for Ai = 1 deg. Note also that, now, a much larger departure from 90 
deg (Ai ~ 0.1 deg) would allow a further reduction of the noise (< 1%). It 
must be noted that for such a high-altitude orbital configuration the gravit- 
omagnetic shift in the out-of-plane direction would amount to 0.10 m yr -1 ; 
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it is certainly a small figure, but, perhaps, not too small if one considers 
that the average shift in the out-of-plane direction of the much lower Mars 
Global Surveyor is 1.6 m after about 5 yr. It poses undoubtedly challenging 
requirements in terms of sensitivity and overall orbit determination accu- 
racy, but it seems not unrealistic to hope that future improvements may 
allow to detect such a small displacement. 

Finally, let us remark that, in principle, also the temporal variations of 
J 2 should be accounted for; however, since at present no secular trends have 
been detected, such changes, mainly seasonal, annual and semi-annual [H] 
would not seriously impact our measurement. 

2.2 The impact of R and GM 

Another martian parameter which must be taken into account is the equa- 
torial radius R along with its uncertainty. Since 



If we assume conservatively for it 5R = 6080 m, i.e. the difference between 
the reference value of MGS95J [11] and the one in p3], the result is depicted 
in Figure [3] for a = 9500 km. If, instead, one takes 5R = 0.04 km the errors 
per degree are as in Figure0]for a = 9500 km. As can be noted, R is the most 
serious source of systematic error, even for high altitudes. Moreover, it must 
also be considered that it would be likely unrealistic to expect improvements 
in the determination of it! by future areocentric missions able to push 5R 
below terrestrial values, i.e. 1 — 0.1 m. Thus, R will likely remain an 
insurmountable obstacle if only one probe is to be used. 

Concerning GM, it turns out that the impact of its mismodelling is of 
no concern being < 1%. 

2.3 Summary 

In summary, in this Section we investigated a nearly polar, high-altitude 
Phobos-like orbit and noted that the Lense-Thirring effect amounts to a 
shift of 0.10 m yr^ 1 in the out-of-plane direction. Concerning the systematic 
errors, the most crucial source of aliasing is the radius R (5fi ~ 10000% 
for 5R = 0.04 km) and, to a lesser extent, the first even zonal harmonic 



Clj oc R l , 



(5) 



the systematic error per degree induced by 5R can be written as 




(6) 
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Figure 3: Systematic percent errors 6fj,j e per degree I due to the uncertainty 
in Mars' radius R, assumed to be 5R = 6080 m, for a = 9500 km, e = 0.01, 
89 deg < i < 91 deg. The errors for £ = 4,6 are as large as 6000% and 
500%, respectively. 
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a= 9500 km, e = 0.01, 6fi = 0.04km 
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Figure 4: Systematic percent errors Sfj,j e per degree I due to the uncertainty 
in Mars' radius R, assumed to be 5R = 0.04 km, for a = 9500 km, e = 0.01, 
89 deg < i < 91 deg. The errors for £ = 4, 6 are as large as 40% and 5%, 
respectively. 
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J2 of the areopotential (5fi < 11%). To reduce them the probe should 
be inserted into an orbit with an inclination close to 90 deg within 10~ 3 
deg or less. In conjunction with such a very tight constraint, it must also 
be hoped that future missions to Mars will improve our knowledge of the 
fundamental parameters of the red planet to a sufficient extent to allow 
for larger departures from the ideal polar geometry: after all, the gravity 
solution MGS95 J represents an improvement of one order of magnitude with 
respect to the previous MGS75D model |15j . However, this may be valid for 
the spherical harmonic coefficients of the areopotential, not for the radius 
of Mars for which an uncertainty of the order of the meter would not be 
enough. We must conclude that the use of only one probe for measuring its 
node Lense-Thirring precession is unfeasible. 



3 Including Phobos and Deimos 

A possible solution to reduce the systematic errors would be to consider a 
linear combination of the nodes of the proposed probe and of the two natural 
satellites of Mars, i.e. Phobos (a = 9380 km, i = 1.075 deg, e = 0.0151) and 
Deimos (a = 23460 km, i = 1.793 deg, e = 0.0002), suitably designed to 
cancel out the impact of just 5J 2 and 5R/R, according to an approach put 
forth for the first time in the context of the terrestrial LAGEOS-LAGEOS 
II test [16]. A possible combination is 

^Deimos + Cl( ^Probe + ^Phobos } ^ 

with 5Q, denoting some sort of Observed-minus-Calculated (O — C) quan- 
tity accounting for any unmodelled/mismodelled features of motion and 
extracted from the data processing with full dynamical force models. By 
purposely leaving the gravitomagnetic force unmodelled, it can be written 
down as 

59, = fiX LT + tl. 2 5J 2 + n R + A : (8) 

here \i = 1 in general relativity, the coefficient Vt ,r is defined as 

^o( obl ) „ . 

n R = R—— = Y / ^iJe, (9) 

1=2 
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and A represents all other remaining mismodelled/unmodelled effects acting 
on the node (e.g. 5J4,5Jq, ...;5GM). By solving for /i one gets 



ci 



('■■2 



^Phobos^Dcimos ^Dcimos^P hobos 

Q Probe oPhob os nPhobosnProbc * 



OD ci mos oProb c rSProbcrSDci mos 

j^ProbcQPhobos^QPhobosQProbe * 



(10) 



Note that c\, contrary to C2, is not defined for i = 90 deg because f2^ robe , 
which are proportional to cosi, vanish for all I at i = 90 deg. The result, 
in principle, would be satisfactory. Indeed, now, the systematic bias of 
J4, Jq, J8-.. on the combination of eq. (J7J) turns out to be <5/i < 6% in 89 deg < 
i < 89.9 deg; the major contribution is due to J4. It is an acceptable result, 
especially, in view of the fact that the admissible range for the inclination 
amounts, now, to about 1 deg and that eq. ([7]) is, by construction, immune 
to the uncertainty in the martian radius; indeed, although it is likely that 
many physical properties and parameters like, e.g., the even zonals, of the 
red planets will be determined with increasing accuracy by the many ongoing 
and planned missions, it is unlikely that the radius will be known to a 
sufficient accuracy to change at an acceptable level the bias induced by it 
(see Figure [3] and Figure H|) . 

The combination of eq. (J7J) is also able to remove the large part of the 
effect of the mismodelling in GM 



1 / 5GM\ 2 " 



2 V GM , 

1=2 



/ D—O 



[ID 



Indeed, the main contribution to eq. (jlip is due to J2, which is canceled out 
by eq. ([7]) with the coefficients of eq. ([TO]) . The other terms in eq. (fTTj) , not 
canceled by the combination of eq. ([7]), yield negligible errors, well below 
1%. 

However, it must be stressed that the node of Phobos undergoes secu- 
lar precessions due to other perturbations of gravitational origin [T7j (non- 
sphericity of Phobos itself, Mars tidal bulge and nutation) which have not 
been considered here and that would affect the combination of eq. ([7]). It 
turns out that the most relevant one is that due to the spherical harmonic 
coefficients C20 and C22 of Phobos itself [18] whose induced secular preces- 
sion amounts nominally to 200 km over 3 yr [T7j . The nutation perturbation 
nominally amounts to 0.3 km over 3 yr [TTJ. The tidally induced precession, 
parameterized in terms of the martian Love number hi would amount nomi- 
nally 0.06 km after 3 vr [17]. but being k 2 = 0.152 ±0.009 pi], its modelling 
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would left a ~ 1 m yr mismodelled trend. The Phobos gravitomagnetic 
shift is 0.1 m yr^ 1 . 

In terms of the detectability of the Lense-Thirring signal with the combi- 
nation of eq. ((TJ), undoubtedly, a major drawback of the strategy of including 
Phobos and Deimos is the poor accuracy with which their orbits can be re- 
constructed with respect to their Lense-Thirring signal. Indeed, while their 
gravitomagnetic shifts are of the order of 1 — 10 cm yr -1 , the latest NASA 
martian ephemerides^l for them yield an accuracy of 1 — 10 km in the radial, 
transverse and out-of-plane orbital components pj_]. Thus, according to the 
present-day level of accuracy including the natural satellites of Mars in the 
combination of eq. (Jj_) would introduce a noise which would overwhelm the 
relativistic trend of interest. 



4 Two dedicated probes 

In principle, the linear combination approach could be followed by using the 
nodes of other two spacecrafts^, although sending to Mars three new probes 
would increase the costs and the difficulties of such a demanding mission. 



4.1 The supplementary orbit configuration 

It is interesting to consider a scenario involving only two probes, named PI 
and P2, in a nearly polar counter-orbiting configuration, proposed for the 
first time by Van Patten and Everitt [20] in the framework of the attempts 
to design a suitable terrestrial mission, and later generalized by Ciufolini [21] 
for other inclinations. By assuming for, say, PI a\ = 9500 km, i\ = 89 deg, 
e\ = 0.01, it turns out that the sum of their nodes would be an observable 
relatively insensitive to departures from the ideal configuration for P2, i.e. 
ci2 = 9500 km, 12 = 91 deg, e2 = 0.01, at least in regard to the mismodelling 
in the even zonals. Indeed, it turns out that for a>2 = 9600 km, 90.05 deg 
< 12 < 92 deg, e2 = 0.03 the bias due to J2 is up to 7%, while the other 
even zonals would have an impact of the order of 1% or less. Unfortunately, 
the uncertainty in the martian radius makes such constraints much more 
stringent. Indeed, S/j, r « 10 5 % for SR = 6080 m and 6/ir « 4000% for 
5R = 0.04 km. Also in this case, the improvements in determining R that 
could realistically be expected from future missions to the red planets may 
be not sufficient. 

4 See on the WEB http://ssd.jpl.nasa.go v/7sat_ephein] 

5 For a polar geometry the Lense-Thirring precession of the pericentre lo of a spacecraft 
vanishes being proportional to cosi. 
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4.2 The linear combination approach 



An approach that could be fruitfully followed within the framework of the 
linear combination strategy with two probes PI and P2, especially in view of 
likely future improvements in the even zonal harmonics of the areopotential, 
consists in designing a combination which, by construction, entirely cancels 
out the bias due to the uncertainty in R 



»Pl 



5R 
~R 



(12) 



being, instead, affected by 5J2, SJ^, 5Jq 
combination of eq. (|7|), with 



It is analogous to the three-nodes 



is 



(13) 



Indeed, as can be noted in Figure [5j 5fij e would be very good (< 1%), 
according to the present-day MGS95J model, provided that inclinations i <C 
90 deg are adopted for the second probe. It turns out that the most relevant 
even zonals would be J4 and Jq with a bias of the order of « 1 — 0.5%. With 
eq. (|13p an inclination range as large as ten degrees would be admissible: it 
is really a great advantage with respect to the other cases examined so far. 



5 Conclusions 

In this paper in examining the possibility of designing a dedicated mission 
to measure the Lense-Thirring effect in the gravitational field of Mars by 
analyzing the secular precessions of the node of one or more spacecrafts we 
focussed on the main systematic errors of gravitational origin and strategies 
to reduce them. 

The main source of systematic bias is the mismodelling in some of the 
parameters entering the multipolar expansion of the Newtonian part of the 
areopotential, especially the Mars' equatorial radius R and, to a lesser but 
non- negligible extent, the even zonal harmonics Ji,£ = 2,4, 6, ... Since the 
resulting aliasing node precessions are proportional to R^a~^ 3 ^ 2+ ^ cos i, a 
high-altitude (a « 9500 km), polar (i = 90 deg) probe would be, in principle, 
an optimal solution, but the consideration of the unavoidable orbit injection 
errors showed that such an option is unfeasible because it would require too 
stringent constraints on the departures from the ideal case i = 90 deg. It 
could be possible, in principle, to suitably combine the nodes of such a probe 
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9600 km, e 2 = 0.03 




Figure 5: Systematic errors Sfxj £ per degree due to the mismodelling 
5J2,SJ4 : ,SJq, ... (MGS95J model) in the uncancelled even zonal harmonics 
in the combination of the nodes of probe PI (a\ = 9500 km, i\ = 89 deg, 
ei = 0.01) and probe P2 (a 2 = 9600 km, 60 deg < i 2 < 89 deg, e 2 = 0.03) 
with eq. (|13p which cancels out entirely the bias due to SR. 
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with those of the natural satellites of Mars, Phobos and Deimos. However, 
the Lense-Thirring shifts of such bodies (1 — 10 cm yr _1 ) are orders of 
magnitude too small with respect to the present accuracy in reconstructing 
their orbits (1 — 10 km). 

A viable option consists in using two probes PI and P2 at high-altitudes 
(di = 9500 km, a2 = 9600 km) and different inclinations {i\ ~ 90 deg, 60 
deg < %2 < 89 deg), and combining their nodes so as to entirely cancel out 
the bias due to R: the resulting bias due to 6J2, 5J4, SJq... would be < 1%, 
according to the present-day MGS95J gravity model, over a range of values 
for the inclination as large as ten degrees. 

A major challenge would certainly be reaching a satisfactory accuracy in 
reconstructing the orbits of such probes whose Lense-Thirring out-of-plane 
shifts would amount to about 10 cm yr _1 ; for example, the spacecraft should 
remain operative for many years, without any failure in the communication 
with the Earth, and it is likely that one or more landers would be required 
as well. 
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